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Abstract. Trajectory hunting (i.e., a technique to find 
air parcels sampled at least twice over the course of a 
few days) is applied to analyze Upper Atmosphere Re- 
search Satellite (UARS) measurements in conjunction 
with the AER photochemical box model. In this study, 
we investigate rapid chlorine activation in the Arctic 
lower stratosphere on 29 Dec 1992 associated with a 
polar stratospheric cloud (PSC) event. Six air parcels 
that have been sampled twice were followed along 5-day 
trajectories at the 465 K (~46 mb) and 585 K (~22 mb) 
levels. A detailed sensitivity study with the AER pho- 
tochemical box model along these trajectories leads to 
the following conclusions for the episode considered: 
(1) model results are in better agreement with UARS 
measurements at these levels if the UKMO temperature 
is decreased by at least 1-2 K; (2) the NAT (nitric acid 
trihydrate) PSC formation scheme produces results in 
better agreement with observations than the STS (su- 
percooled ternary solution) scheme; (3) the model can 
explain the UARS measurements at 585 K, but under- 
estimates the CIO abundance at 465 K, suggesting some 
inconsistency between the UARS measurements at this 
level. 


Introduction 

To understand the current and predict the future ozone 
depletion in the polar regions it is necessary to separate 
the effects due to dynamical and chemical processes. 
One effective strategy for achieving this goal is to iden- 
tify air parcels that have been sampled twice over the 
course of a few days. The difference between initial 
and final concentrations of ozone (or any other species) 
for these parcels can be attributed solely to chemical 
processing since dynamical effects are eliminated. This 
Lagrangian approach has been used by von der Gathen 
et ai [1995] and Rex et aL [1998] to quantify chemical 
ozone depletion from ozonesondes launched at differ- 
ent points along a parcel trajectory (“matches”) in the 
Arctic. Becker et ai [1998] used a photochemical box 
model to analyze the measured O3 loss along matched 
parcels. This technique can also be applied to mea- 
surements obtained by different platforms. For exam- 
ple, Pierce et ai [1997] compared trajectory-matched 
observations from the UARS Halogen Occultation Ex- 
periment (HALOE) with aircraft ER-2 measurements. 
Bacmeister et ai [1998] applied trajectory mapping 
to a comparison of long-lived tracer data obtained by 
the ER-2 and CRyogenic Infrared Spectrometers and 
Telescopes for the Atmosphere (CRISTA) in November 
1994. 

Analysis of short-lived species along parcel trajectories 
requires a photochemical model. Utilization of UARS 
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Cryogenic Limb Array Etalon Spectrometer (CLAES) 
and Microwave Limb Sounder (MLS) data combined 
with a photochemical box model looks particularly 
promising in this context due to a relatively large num- 
ber of measured species (0 3 , HNO3, CIO, CIONO2, 
N 2 0, CH 4 , H 2 0) and very frequent sampling (~1300 
profiles per day). The goal of this study is to compare 
the AER box model representation of photochemical 
and PSC processes with the UARS measurements for 
air parcels sampled twice and to check the consistency 
between CLAES and MLS data. 

Identification of matched trajectory 

We examine an episode of rapid chlorine activation 
associated with PSC activity in the Arctic polar re- 
gion. We concentrated on events for which the evi- 
dence of chlorine activation is unambiguous, i.e., day- 
time MLS measurements showing unenhanced values of 
CIO immediately prior to formation of a PSC. In one 
such episode, the CIO concentration increased from near 
zero on 24 Dec 1992 to 1.4 (0.8) ppbv on 29 Dec 1992 
at 465 K (585 K) (see Figure 1). For this case, 400 
equally spaced parcels were initialized in a 20° x 40° 
latitude/longitude box encompassing the enhanced CIO 
points on 29 Dec and 5-day isentropic back trajecto- 
ries were calculated using the model described by Man- 
ney et al. [1994] (with horizontal winds and tempera- 
tures from the UK MO data assimilation). Because di- 
abatic effects are negligible over a 5-day interval in the 
lower stratosphere, cross-isentropic flow was ignored. 
From this original set of parcels, we identified candi- 
date parcels that were nearly coincident (in space and 
time) with enhanced CIO points on 29 Dec and with 
any MLS measurement points 5 days earlier on Dec 24. 
Also, both initial and final points of these parcels were 
within the vortex (based on UKMO potential vortic- 
ity gradients). The number of matches depends on the 
number of individual enhanced CIO points, the num- 
ber of parcels initialized in the original box, the me- 
teorological conditions, and of course the “match” cri- 
teria. For example, sensitivity studies show that for 
the episode considered, the number of matches drops 
from 52 (35) at 465 K (585 K) for At < 2 h, Alatitude 
< 4°, Alongitude < 350 km to 14 (6) for At < 2 h, 
Alatitude < 2°, Alongitude < 200 km. Many matches 
are duplicative in the sense that several parcels may 
start off near the same high-ClO point, follow similar 
paths, and thus also end up nearest the same point at 
the termination of the back trajectory calculation. For 
the results discussed below, we have selected a subset 
of 6 parcels on each level that represents the widest 
range of behavior possible; ail selected parcels satisfy 
the stricter match criteria (At < 2 h, Alatitude < 2°, 
Alongitude < 200 km). Other uncertainties associated 
with the match technique are discussed in Rex et al. 
[1998]. Note that this technique should only be used 
for simulations shorter than ~10 days because of the 
large uncertainties in trajectory calculations done for 
longer periods. 



Model Initialization 


The AER photochemical box model [Danilin et al. , 
1998] with NAT and STS PSC schemes and the recom- 
mended list of heterogeneous reactions [DeMore et al. , 
1997] is used. The NAT surface area density (SAD) is 
calculated assuming a particle diameter of 1 /mi with 
a density of 1.62 g/cm 3 , while the STS SAD is defined 
according to [ Tabazadeh et al., 1994] keeping the par- 
ticle concentration constant. The values of the MLS 
and CLAES species are obtained at the exact initial 
and final parcel positions through interpolation. The 
model is initialized using MLS version 4 (v.4) O3 (pro- 
files and column above 100 hPa) and HNO3, prototype 
non-linear retrievals of MLS H 2 0 [Pumphrey, 1998], and 
CLAES version 8 (v.8) C10N0 2 , N 2 0, and CH4. To- 
tal inorganic chlorine, Cl y , is defined from the Cl y -N20 
correlation according to W oodbridge et al. [1995] using 
CLAES measurements of N 2 0 at the initial points; this 
approach produces 2. 0-2.5 ppbv and ~2.9 ppbv of CL at 
465 K and 585 K, respectively. For all trajectories, CIO 
is initialized to C10 MLS . Since no HALOE HC1 mea- 
surements are available for the points of interest, initial 
HCl is set to Cl y - C10N0!? LAES - C10 MLS . Total in- 
organic bromine, Br y , is defined according to Wamsley 
et al. [1998] with initial Br0=Br0N0 2 =Br y /2. Since 
no reliable NO z (= N0+N0 2 ) measurements are avail- 
able for the points considered, we perform a sensitivity 
study in which NO x is set to 2 pptv or 0.2 ppbv (based 
on aircraft measurements [Kawa et al. , 1992]), and we 
assume that initial N0=N0 2 =N0 z /2. Total reactive 
nitrogen, NO„, is taken to be HNO^ LS + C10N0£ LAES 
+ NO z . Sulfate aerosol surface area is initialized ac- 
cording to Massie et al. [1998] using the CLAES aerosol 
extinction measurements at 780 cm -1 . 

Results of trajectory calculations 

MLS detected a rapid increase in CIO of 0.5-1 ppbv in 
the lower stratosphere above northern Canada between 
24 Dec and 29 Dec 1992 (Figure 1). Six air parcels 
that satisfy the match criteria described above are also 
shown in Figure 1, together with the temperature and 
solar zenith angle (SZA) along their trajectories. Tem- 
peratures dipped below 195 K in the 1-2 days preced- 
ing the MLS observation of high CIO, suggesting that 
PSC processing converted chlorine into Cl 2 and C1 2 0 2 , 
which were photolyzed to CIO at the final points (where 
SZA<90°). 

Figure 2 compares model results and data for CIO and 
C10N0 2 along the trajectory of parcel #A at 585 K 
and parcel #a at 465 K. The blue curves depict the 
behavior for the NAT PSC scheme and UKMO temper- 
ature. It is obvious that this approach cannot explain 
the CIO measurements. Manney et al. [1996] showed 
that, compared to radiosondes, UKMO temperatures 
have a systematic warm bias of ~l-3 K during most of 
the northern winter. To explore the possible influence 
of such a temperature bias, we show model results for 
the NAT scheme with the UKMO temperature lowered 



by 2 K. In this case (grey lines), the model agrees with 
the CIO measurements at 585 K, but overestimates the 
CIONO2 depletion. By additionally increasing the ini- 
tial C10N0-> by 15% [ Mergenthaler et ai, 1996] (while 
keeping Cl y constant), the model (green lines) can be 
made to agree with the measurements at 585 K. How- 
ever, the model still underestimates CIO and CIONO2 
at 465 K even with all these adjustments. For the STS 
scheme (red lines), the model underestimates CIO at 
both levels. The CLAES v.8 aerosol extinction mea- 
surements (converted to SAD) are consistent with the 
model calculations in the final points. However, these 
measurements are not useful for distinguishing between 
the NAT ans STS schemes during the episode consid- 
ered since the PSC event was not captured in the fi- 
nal points of the matched trajectories. Very similar 
results are obtained for the other parcels and are not 
shown here. Comparisons of model results and UARS 
measurements are not shown for O3 and HNO3 since 
the differences between their initial and final values are 
small and generally consistent with the model within 
the range of the measurement uncertainties. 

Table 1 summarizes our model sensitivity analysis for 
all parcels. The model can be made to agree with the 
CIO and CIONO2 measurements along each parcel tra- 
jectory at 585 K by varying the temperature or the ini- 
tial CIONO2 or NO x within reasonable ranges. How- 
ever, there is a persistent model underestimate of CIO 
at 465 K for all considered sensitivity runs. Previous 
model attempts also underestimated the MLS CIO mea- 
surements at 465 K [Chipperfield et ai, 1996; Lutman 
et al. } 1997]. 

Discussion 

The failure of the straightforward and well-constrained 
comparison in our case raises the question of whether 
the MLS v.4 CIO and CLAES v.8 CIONO2 measure- 
ments are consistent at 465 K during/after PSC events. 

To answer this question, we estimate the maximum 
CIO concentrations allowed by stoichiometry arguments 
based on the following assumptions: (1) the CIONO2 
+ HC1 — » CI2 -F HNO3 reaction is a principal mech- 
anism of chlorine activation, (2) initial [HC1] > initial 
[CIONO2], (3) initial [Cl 2 ] and [C1 2 0 2 ] are very small, 
and (4) all active chlorine is in CIO at the final point. 
The first assumption is very reasonable, bearing in mind 
the temperature history along the trajectories and the 
temperature dependence of the C10N0 2 + H 2 0 and 
HOC1 + HC1 reactions [DeMore et ai , 1997], The sec- 
ond assumption is consistent with N 2 OCLAES-derived 
Cly for the episode studied. The third assumption is 
supported by the warm temperature along the back- 
ward trajectories and the CIO and C10N0 2 measure- 
ments 10 days prior to 24 Dec 1992, suggesting that 
conversion of HC1 and C10N0 2 to other active chlorine 
forms (C1 2 0 2 and Cl 2 ), undetectable by MLS, is un- 
likely. The fourth assumption is too strong for the final 
points with SZA of 85°-90°, but it provides an upper 
bound for CIO. Based on these assumptions, the max- 



irnum increase in chlorine monoxide, AClO max , equals 
2xC10N0 2 im “ a< - C10N0^' iai (since the C10N0 2 + 
HC1 reaction produces CI 2 until a complete consump- 
tion of the smaller reactant (i.e. C10N0 2 mtt * a/ ): and 
since a part of activated chlorine is converted back 
to CIONO^ 1 ™*). Taking into account the 15% preci- 
sion of the CLAES CIONO? measurements, one gets 

AC10 M/lx = 2.3xC10N0j ni ‘ ial - 0.85xClONOp ,a '. 
A possibility that the whole initial HC1 amount could 
be consumed requires both prolonged PSC processing 
in a presence of sunlight and high NO x (several ppbv). 
However, this possibility can be ruled out based on the 
SZA history along these parcels (see Figure 1) and rel- 
atively low NO x (up to 0.2 ppbv) [Kawa et al. , 1992]. 

Table 2 compares the changes in CIO according to 
the MLS measurements and our stoichiometric calcu- 
lations. The MLS v.4 CIO and CLAES v.8 C10N0 2 
data are consistent within the range of their uncertain- 
ties if AC10 MLS -2or < AC\O max . Table 2 shows that 
these measurements are inconsistent for parcel #b. It 
is unlikely that these measurements are consistent for 
parcel #e (where about half of the active chlorine is in 
the CI 2 O 2 form at the final point, for which SZA^89°). 
Analysis of a larger subset of UARS data is required 
to elucidate whether MLS V.4 overestimates CIO or 
CLAES v.8 overestimates C10N0 2 after PSC events, 
or whether the model fails to adequately simulate other 
physical processes. Various refinements to the MLS and 
CLAES retrieval algorithms are ongoing and future ver- 
sions of these data sets may exhibit improved consis- 
tency with model calculations. 
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Figure 1. MLS v.4 CIO (ppbv) at 465 K (left) and 
585 K (right) on 29 Dec (top row) and 24 Dec 1992 
(second row). Isentropic 5-day back trajectories for 6 
air parcels that were sampled twice are overlaid on the 
24 Dec plots in the second row, with thick dots cor- 
responding to the locations of the MLS measurements 
above Canada on 29 Dec that indicate enhanced CIO 
(defined as >1.0 ppbv at 465 K and >0.75 ppbv at 
585 K). The evolution of temperature (third row) and 
SZA (bottom row) along these trajectories are shown 
using the same colors as for the parcels. 


Figure 1 . MLS v.4 CIO (ppbv) at 465 K (left) and 585 K (right) on 29 Dec (top row) and 24 Dec 1992 (second 
row). Isentropic 5-day back trajectories for 6 air parcels that were sampled twice are overlaid on the 24 Dec plots 
in the second row, with thick dots corresponding to the locations of the MLS measurements above Canada on 
29 Dec that indicate enhanced CIO (defined as >1.0 ppbv at 465 K and >0.75 ppbv at 585 K). The evolution of 
temperature (third row) and SZA (bottom row) along these trajectories are shown using the same colors as for 
the parcels. 


Figure 2. Comparison of MLS v.4 CIO and CLAES 
v.8 CIONO 2 and aerosol (converted from the 780 cm -1 
extinction to surface area density, SAD) (black sym- 
bols) with AER model calculations along the trajec- 
tory of parcels #A and #a (denoted by red lines in 
Figure 1) at 585 K (top) and 465 K (bottom), respec- 
tively. Blue lines: NAT scheme, UKMO temperature; 
grey lines: NAT scheme, UKMO temperature lowered 
by 2 K; green and red lines: UKMO temperature low- 
ered by 2 K, initial CIONO 2 set to l.lSxCiONO^ LAES j 
NAT and STS schemes, respectively. The grey and 
green lines coincide in the SAD plot. 


Figure 2. Comparison of MLS v.4' CIO and CLAES v.8 C10N0 2 and aerosol (converted from the 780 cm -1 
extinction to surface area density, SAD) (black symbols) with AER model calculations along the trajectory of 
parcels #A and #a (denoted by red lines in Figure 1) at 585 K (top) and 465 K (bottom), respectively. Blue 
lines: NAT scheme, UKMO temperature; grey lines: NAT scheme, UKMO temperature lowered by 2 K: green 
and red lines: UKMO temperature lowered by 2 K, initial CIONO 2 set to 1.15xC10N02 :LAES , NAT and STS 
schemes, respectively. The grey and green lines coincide in the SAD plot. 
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Table 1 . Comparison of model calculations of CIO and C10N0 2 at the final points with UARS measurements. 
“T” and “T-2K” denote the NAT scheme with T UKMO and T UKM D - 2 K, respectively; “STS” denotes the STS 
scheme with T UKMO - 2 K; “fit” denotes the NAT scheme with T UKMO - 2 K and initial CIONO 2 and NO x 
varying from CIONO? IjAI * :S to l.lbxClONO^ LAES and from 0.002 to 0.2 ppbv, respectively. 


Species 

Par- 

cel 

465 K 

UARS, 

ppbv 

T, 

NAT 

Model 

T-2K, 

NAT 

T-2K, 

STS 

Par- 

cel 

585 K 

UARS, 

ppbv 

T, 

NAT 

Model 

T-2K, 

NAT 

T-2K, 

STS 

fit, 

NAT 

CIO 

a 

1.10T0.48 

0.07 

0.44 

0.18 

A 

0.79±0.38 

0.24 

0.78 

0.24 

0.41 

CIONO 2 


0.64±0.10 

0.59 

0.36 

0.51 


0.94±0.14 

0.88 

0.68 

0.88 

0.88 

CIO 

b 

1.28T0.49 

0.39 

0.58 

0.42 

B 

0.82±0.38 

0.20 

0.24 

0.20 

0.43 

C10N0 2 


0.33T0.05 

0.30 

0.04 

0.15 


0.87±0.13 

0.90 

0.91 

0.92 

0.86 

CIO 

c 

1.39±0.49 

1.02 

1.05 

1.02 

C 

0.88±0.39 

0.21 

0.26 

0.21 

0.51 

C10N02 


0.50T0.08 

0.07 

0.01 

0.05 


0.95±0.14 

0.94 

0.95 

0.95 

0.87 

CIO 

d 

1.39±0.49 

0.86 

0.97 

0.72 

D 

0.81±0.38 

0.26 

1.26 

0.23 

0.51 

C10N02 


0.65T0.10 

0.17 

0.02 

0.28 


1.21±0.18 

0.86 

0.40 

0.88 

1.17 

CIO 

e 

1.49±0.49 

0.71 

0.71 

0.71 

E 

0.78±0.38 

0.19 

0.55 

0.20 

0.57 

C10N02 


0.37T0.06 

0.02 

0.01 

0.00 


1.18±0.18 

0.95 

0.81 

0.95 

1.00 

CIO 

f 

1.49±0.49 

0.79 

0.82 

0.76 

F 

0.78±0.38 

0.09 

0.28 

0.09 

0.41 

C10N02 


0.40±0.06 

0.07 

0.01 

0.01 


1.10±0.17 

0.96 

0.90 

0.96 

0.97 


Table 2. Change in CIO (ppbv) between final and 
initial points at 465 I< according to MLS v.4 data and 
stoichiometric calculations (see text). 


parcel 


MLS 

Stoichiometry 


ACIO 

AC10-2fr 

AC10 mftI 

&C\0 MAX 

a 

1.11 

0.54 

0.58 

0.86 

b 

1.28 

0.71 

0.47 

0.64 

c 

0.82 

0.25 

0.68 

0.94 

d 

0.82 

0.25 

0.49 

0.76 

e 

0.92 

0.35 

0.39 

0.56 

f 

0.92 

0.35 

0.72 

0.94 
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Figure 1 . MLS v.4 CIO measurements (in ppbv) at 465 K (left column) and 585 K (right column) 
on Dec 29 (top row) and Dec 24, 1992 (second row). Six air parcels sampled twice are shown in the 
second row with thick dots (final points) corresponding to the locations of MLS measurements on 
Dec 29, 1992 showing high CIO content above Canada (red dots at 465 K and yellow dots at 585 K). 
Evolutions of temperature (third row) and SZA (low row) along these parcels are also shown by the 
same colors as the parcels. 
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Figure 2. Comparison of MLS v.4 CIO and CLAES v.8 CIONO 2 and aerosol (converted 
from the 780 cm -1 extinction to SAD) measurements (black symbols) at 465 K (bootom 
line) and 585 K (top line) with AER model calculations for the parcels #a and #A (i.e. 
shown by red lines in Figure 1). Blue lines: NAT scheme, UKMO temperature; grey 
lines: NAT scheme, UKMO temperature lowered by 2 K; green and red lines: UKMO 
temperature lowered by 2 K, initial CIONO 2 equals to 1.15C10N02 ae,s , NAT and STS 
schemes, respectively. The grey and green lines coincide at the SAD plot. 









